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We describe the temperature-chemical potential phase diagrams of holographic models of a range 
of strongly coupled gauge theories that display chiral symmetry breaking/restoration transitions. 
The models are based on the D3/probe-D7 system but with a phenomenologically chosen running 
coupling/dilaton profile. We realize chiral phase transitions with either temperature or density that 
are first or second order by changing the dilaton profile. Although the models are only caricatures 
of QCD they show that holographic models can capture many aspects of the QCD phase diagram 
and hint at the dependence on the running coupling. 



I. INTRODUCTION 

The QCD phase diagram is notoriously difficult to 
compute. Firstly the physics associated with deconfine- 
ment or chiral symmetry restoration is strongly coupled 
where we traditionally do not know how to compute. Sec- 
ondly at finite density lattice gauge theory, the first prin- 
ciples simulation of the theory on supercomputers, suffers 
a "sign problem" that means Monte Carlo methods break 
down. In fact with light quarks there is no clear order pa- 
rameter for deconfincment so we will concentrate on the 
chiral transition. Progress has been made by identifying 
effective theories of the transitions and through lattice 
computations at low density. [1] provides a review of 
the standard picture. It is believed for QCD, with the 
physical quark masses, that the phase transition with 
temperature is a smooth cross over (becoming a second 
order transition as the up and down quark masses go to 
zero). At zero temperature the transition with density 
is believed to be first order. There must therefore be a 
critical point where the first order line ends in the tem- 
perature density plane. 

In the last ten years the AdS/CFT Correspondence [2- 
4] and more general gauge/gravity dualities have emerged 
as a new tool for the study of strongly coupled problems. 
It is interesting to ask whether these holographic models 
can in principle describe a phase diagram like that of real 
QCD. One recent attempt in this direction can be found 
in [5] where a 5d holographic model of a strongly coupled 
gauge theory with a running coupling was shown to give 
an appropriate looking phase diagram. In that model 
the order parameter is associated with confinement and 
it does not include quark fields. In true QCD the order 
parameter across the phase diagram is the quark bilinear 
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condensate. In this paper we will study a holographic 
model of quark fields and again attempt to reproduce 
key features of the QCD phase diagram. 

Quarks can be introduced into the AdS/CFT Corre- 
spondence through probe branes [6-10] and a number of 
systems with chiral symmetry breaking have been devel- 
oped [11-16]. We can not of course describe true QCD 
holographically because the dual, if it exists, is not known 
and is probably very complicated (and strongly coupled, 
at least, in the UV). Our analysis is therefore in the spirit 
of AdS/QCD [17, 18], a phenomenological modeling of 
the QCD phase diagram. If one could model the phase 
diagram correctly one might hope to then predict other 
features of the theory such as time dependent dynamics 
during transitions and so forth. 

Our models will be in the context of the simplest brane 
construction of a 3+ld gauge theory with quarks which is 
the D3 /D7 system of Fig 1 [6-10] . The basic gauge theory 
is large N, Af = A super Yang-Mills with Nf quark fields. 
We will work in the quenched approximation where we 
neglect quark loops. On the gauge theory side we do 
not backreact the D7 branes, that provide the quarks, on 
the geometry but instead work in the probe approxima- 
tion [8]. The theory has a.U{l) symmetry under which a 
fcrmionic quark anti-quark condensate has charge 2 and 
plays the role of C/(l) axial [11]. Although the theory 
does not have a non-abelian chiral symmetry (as for ex- 
ample the Sakai Sugimoto model does [16]) this is not 
important in the quenched approximation since the dy- 
namics of the formation of the quark condensate is flavour 
independent. The model is very simple to work with hav- 
ing a background geometry that is just AdSs x and 
the gauge theory is 3-|-ld at all energy scales. 

So called top down models of this type exist with chiral 
symmetry breaking. Supergravity solutions exist that 
correspond to the AdS space being deformed in reaction 
to a running coupling introduced by a non-trivial dilaton 
profile [11, 13]. In cases where the coupling grows in the 
infra-red (IR), breaking the conformal symmetry, chiral 
symmetry breaking is induced. These models have very 




FIG. 1: A schematic of the D3/D7 showing our conven- 
tions. The D3-D3 strings generate the A/" = 4 theory, the 
D3-D7 string represent the quarks and D7-D7 strings describe 
mesonic operators. 



specific forms for the running coupling and are typically 
singular somewhere in the interior. At the string theory 
level a full interpretation is lacking. 

A yet simpler and completely computable case with 
chiral symmetry breaking is provided by introducing a 
background magnetic field associated with U{1) baryon 
number [15]. Such a background source can be described 
by a gauge field on the surface of the D7 brane. A chiral 
condensate is induced. Very simplistically one can think 
of the B-field as introducing a scale that breaks the con- 
formal symmetry as the strong coupling scale Aqcd does 
in QCD allowing the strong dynamics to form the quark 
condensate. 

In recent papers [19-23] we have explored the phase 
structure of the theory with magnetic field. Tempera- 
ture can be introduced through an AdS Schwarzschild 
black hole in the geometry [4] . Density and chemical po- 
tential can be added through the temporal component of 
the U{1) baryon number gauge field [25-27]. The phase 
diagram [19] is shown in Fig 2(a). The chiral restoration 
transition was found to be first order with temperature 
and second order with density. A critical point lies be- 
tween these regimes. In addition there is an extra tran- 
sition associated with the formation of density and the 
mesons of the system melting into the background plasma 
[28, 29]. In places in phase space this merges with the 
chiral transition but in other places it separates and can 
be either first or second order. Such a phase with a quark 
density but chiral symmetry breaking could potentially 
exist in QCD. It would be nice to also describe this tran- 
sition as a deconfincment transition but firstly the A/" = 4 
background does not induce linear confinement and sec- 
ondly the presence of any temperature leads to screening 
of the quarks at the length scale of the inverse temper- 
ature. The meson bound states are closer in spirit to 
atomic bound states than QCD-like mesons. Neverthe- 
less they are being disrupted by the background plasma 
so the existence of a phase with melted mesons but chiral 
symmetry breaking at least leads one to speculate on a 
possible separation of deconfinement and chiral restora- 
tion behaviour in the QCD phase plane. 

In the gravity picture of the transitions the crucial 
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{a)T-ii phase diagram. The inset diagrams are 
representative probe brane embeddings (dotted lines), 
where a black disk represents a black hole. 
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FIG. 2: The phase diagrams of the massless J\f — 2 gauge 
theory with a magnetic field. First order transitions are shown 
in blue, second order transitions in red. The temperature is 
controlled by the parameter T, chemical potential by /i and 
electric field by E. 
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FIG. 3: The phase diagram of the massless axion/dilaton 
gauge theory in [31]. 



transitions are between the three sorts of D7 brane em- 
bedding shown in Fig 2(a). An embedding that curves 
in the holographic space to miss the black hole repre- 
sents the chiral symmetry breaking phase. These solu- 
tions have stable and discrete linearized fluctuations that 
correspond to the mesons of the theory. The second type 
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of embedding curves off axis but ends on the black hole. 
It describes a phase with chiral symmetry breaking but 
the linearized fluctuations are now replaced by in-falling 
quasi-normal modes of the black hole which describe un- 
stable mesonic fluctuations of a quark plasma. Finally 
the flat embedding ending on the black hole describes 
melted mesons and chiral symmetry restoration. The 
transitions occur when the magnitude of the D7 action 
for these different cases cross. 

How generic to the holographic description is this 
phase structure? Keeping within the top down analy- 
sis one can change parameters and sec what effect they 
have on the phase diagram. For example in [21] we traded 
density for an electric field parallel to the magnetic field. 
The electric field tries to dissociate mesons by accelerat- 
ing the quark and anti-quark in opposite directions and 
so opposes the formation of a chiral condensate. The 
phase diagram of the theory in the T-E plane at fixed B 
is also shown in Fig 2(b). The chiral transitions are again 
first order with temperature but second order with den- 
sity although there is some change in the meson melting 
transition order^. We followed up this analysis in [22] to 
consider the case of mutually perpendicular electric and 
magnetic fields - there the chiral transition is first order 
in nature throughout the full T-E plane. 

A recent paper [31] performed a similar analysis with 
a running dilaton geometry. The geometry is that of [32] 
in which there is a non-zero profile for both the dilaton 
and axion fields in AdS. The field theory is the A/" = 4 
gauge theory with a vev for both Tri^^ and TrFF which 
preserves supersymmetry at zero temperature but dis- 
plays confinement. A D7 was introduced in a supersym- 
metry breaking fashion and chiral symmetry breaking is 
observed. The temperature density chiral transformation 
was first order throughout the plane and is shown in Fig 
3. It shows the same three phases as the magnetic field 
case. An extra component of the analysis in [31] was to 
note that in the confining geometry with a running dila- 
ton a baryonic phase was also present. A baryon vertex 
is described by a D5 brane wrapped on the of the 
AdSs X space. In the pure Af = 4 theory such vertices 
shrink to zero. However in the running dilaton geometry 
the large IR value of the dilaton stabilizes the D5 embed- 
ding. Solutions exist that link the D5 to the D7 brane 
embedding with a balancing force condition. These con- 
figurations describe the gauge theory with finite baryon 
density rather than finite quark density. This phase sets 
in at finite chemical potential and then persists to infinite 
chemical potential (as shown in Fig 3) which is certainly 
unlike QCD. We will not focus on this phase in this paper 
but it would be interesting to study it in future work to 
find models that have a baryonic phase in some interme- 
diate range of chemical potential like QCD. 



^ There may be an instability in Fig 2 due to the WZ term contri- 
bution [22, 30]. 



These phase structures are very interesting and sur- 
prisingly complex but do not match the expectations in 
QCD. In QCD we need a second order transition with 
temperature and a first order transition with density to 
the chirally symmetric phase. 

Here we want to work in a much more generic frame- 
work to ask what phase structures it is possible to get 
in the holographic description and to try to force our- 
selves onto a representation of the QCD phase diagram. 
We will therefore take a bottom up approach within the 
model and allow ourselves to dial the running of the 
gauge coupling by hand. We will have a dilaton profile 
that smoothly transitions from a UV conformal regime 
to an IR conformal regime through a step of variable 
height and width. Such an ansatz allows one to consider 
runnings that range from precocious growth in the IR 
to more walking like dynamics [33]. We used a similar 
ansatz in previous papers to study the impact of walking 
[14] on meson physics and as a mechanism for generating 
infiation [34]. Here, in a completely new analysis of the 
phase structure of these models, we find that with the 
simple step ansatz we can move from a totally first order 
transition in the phase plane to a configuration similar to 
that we obtained with a B field (a first order transition 
with temperature but second order with density). With 
this ansatz we can not achieve a second order transition 
with temperature. 

The model directly suggests other phenomenological 
generalizations though. In particular, if we think of the 
running dilaton profile as a short cut for including the 
back reaction due to the quark fields/D7 brane ^, then 
it is natural to break the 5*0(6) symmetry of AdSs in 
the dilaton in the same fashion as the D3/D7 system's 
geometry. This allows us an extra phenomenological free- 
dom to distort the dilaton or black hole horizon. These 
simple changes do allow us to reproduce a wide range of 
phase diagrams including QCD-likc ones as we will show 
below. We will discuss the simple geometric reasons for 
the emergence of first or second order transitions in these 
different scenarios. 

Our conclusion therefore is that the holographic model 
has no intrinsic problem with mimicking the QCD phase 
diagram and these systems may therefore be phcnomcno- 
logically useful in the future. 

II. THE MODEL 

First let us review the gravity dual description of the 
symmetry breaking behaviour of our strongly coupled 
gauge theory. 



^ In [35], D3-D7 solutions at finite temperature and chemical po- 
tential, with the inclusion of dynamical flavor effects, have been 
derived and studied as full-fledged top-down models. They will 
be the first step towards the top-down study of phase transitions 
in D3-D7 systems with dynamical flavors. 
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Dp-branes are p dimensional membrane like objects 
to which the ends of open strings are tied. The weak 
coupling picture for our D3/D7 set up is shown in Fig 
1 [6-10] - there are TV D3 branes and the lightest string 
states with both ends on the D3 generate the adjoint 
representation fields of the A/" = 4 gauge theory. Strings 
stretched between the D3 and the D7 arc the quark fields 
lying in the fundamental representation of the SU {N) 
group (they have just one end on the D3). 

In the strong coupling limit the D3 branes in this pic- 
ture are replaced by the geometry that they induce. We 
will consider a gauge theory with a holographic dual de- 
scribed by the Einstein frame geometry AdSs x 



where Ty = (27r) 



^QuY and Ty = when we 



^2 jj)2 



g'^dnl + dwl + dwl) , (1) 



where we have split the coordinates into the xa+i of the 
gauge theory, the g and which will be on the D7 brane 
world-volume and two directions transverse to the D7, 
W5,wq. The radial coordinate, = g'^ + w'^ + Wq, cor- 
responds to the energy scale of the gauge theory. The 
radius of curvature is given by i?4 = AirgsNa"^ with N 
the number of colours. The r — > oo limit of this the- 
ory is dual to the A/" = 4 super Yang-Mills theory where 
9s = 5uv constant large r asymptotic value of the 

gauge coupling. 

In addition we will allow ourselves to choose the profile 
of the dilaton as r — )■ to represent the running of the 
gauge theory coupling, e'^ = /3, where the function /3 1 
as r — >■ oo. An interesting phenomenological case is to 
consider a gauge coupling running with a step of the form 
[14, 34] 



^ + 1 - Atanh[r(r - A)] 



(2) 



Of course in this case the geometry is not back reacted 
to the dilaton and the model is a phenomenological one 
in the spirit of AdS/QCD[17, 18]. This form introduces 
conformal symmetry breaking at the scale A = A/27ra' 
which triggers chiral symmetry breaking. The parameter 
A determines the increase in the coupling across the step. 

We will introduce a single D7 brane probe [8] into the 
geometry to include quarks - by treating the D7 as a 
probe we are working in a quenched approximation al- 
though we can reintroduce some aspects of quark loops 
through the running coupling's form if we wish (or know 
how). This system has a U{1) axial symmetry on the 
quarks, corresponding to rotations in the w^-wq plane, 
which will be broken by the formation of a quark con- 
densate. 

In the true vacuum at T = the brane will be 
static. We must find the D7 embedding function e.g. 
W5{g),we = 0. The Dirac Born Infold (DBI) action in 
Einstein frame is given by 



Sd7 = -T7 J d''Ce'^^-detP[GU 



^-Tj / dV dg g'pJl + [d,W5)\ 



(3) 



have integrated over the 3-sphere on the D7. The equa- 
tion of motion for the embedding function is therefore 



l3g^dgW5 



^2w,g^^l + {d,w,Y^ = Q. (4) 



The UV asymptotic of this equation, provided the dilaton 
returns to a constant so the UV dual is the A/" = 4 super 
Yang-Mills theory, has solutions of the form 



g^ 



(5) 



where we can interpret m as the quark mass (ruq = 
m/2TTa') and c is proportional to the quark condensate. 

The embedding equation (4) clearly has regular solu- 
tions W5 = m when gyM independent of r - the flat 
embeddings of the Af = 2 Karch-Katz theory [8] . Equally 
clearly if dp/dr^ is none trivial in then the second 
term in (4) will not vanish for a flat embedding. 

There is always a solution wc, = Q which corresponds 
to a massless quark with zero quark condensate (c = 0). 
In the pure Af — 2 gauge theory with (3 — 1 this is the 
true vacuum. In the symmetry breaking geometries this 
configuration is a local maximum of the potential. 

If the coupling is larger near the origin then the D7 
brane will be repelled from the origin ^ . The parameter 
r spreads the increase in the coupling over a region in r 
of order T^^ in size. 

We display the embeddings for some particular cases 
in Fig 4. Note that we have chosen parameters here that 
make the vacuum energy of the theory the same in each 
case. The vacuum energy is given by minus the DBI 
action evaluated on the solution. In fact this energy is 
formally divergent corresponding to the usual cosmolog- 
ical constant problem in field theory. As usual we will 
subtract the UV component of the energy to renormal- 
ize. 

The symmetry breaking of these solutions is visible 
directly [11]. The U{1) symmetry corresponds to rota- 
tions of the solution in the w^-wq plane. An embedding 
along the g axis corresponds to a massless quark with 
the symmetry unbroken (this is the configuration that is 
preferred at high temperature and it has zero condensate 
c). The symmetry breaking configurations though map 
onto the flat case at large g (the UV of the theory) but 
bend off axis breaking the symmetry in the IR. 

One can interpret the D7 embedding function as the 
dynamical self energy of the quark, similar to that emerg- 
ing from a gap equation [14]. The separation of the D7 



^ In fact there is a competition between the increased action from 
the D7 entering the region with larger dilaton and the derivative 
cost of the D7 bending to avoid it. This leads to a critical value of 
A to trigger chiral symmetry breaking. For example for A = 1.7 
and r = 1 Ac = 2.1. In this paper we will consider only super- 
critical values of A. 
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FIG. 4: Example coupling flows (2) (top) and the induced 
D7 brane embeddings/quark self energy (bottom) with the 
parameter choices shown in the table. 



hole is the natural eandidate since it has intrinsic ther- 
modynamic properties such as entropy and temperature. 
We make the coordinate transformation [11] 



rdr 



(r4 - r|^)i/2 



2w^ = 



dw 
w 



with \/2wH = rn, such that the metric becomes 

„,,2 



i?2 



(dp2 + p^dnl + AL^ + L^dn\), 



where 
9t 



(w^ — w|f )^ 



gx 



4 I 4 
W + W 



H 



w = \J + L"^ , p ~ w sin 6 , L ^ w cos ( 



Now we have to transform /? also 
and therefore 



4 I 4 
W + Wjj 



fi = A + l- A tanh 



- A 



(8) 



(9) 



(10) 



(11) 



(12) 

Note that for wh — >■ 0: w ^ r, p ^ g and L ^ w^, if we 
set wq = 0. 



from the g axis is the mass at some particular energy 
scale given by g. 



A. Temperature 

Temperature can be included in the theory by using the 
AdS-Schwarzschild black hole metric. In Einstein frame 
we have 



B. Chemical Potential 

The DBI action for the D7 brane naturally includes 
a surface gauge field which holographically describes the 
quark bilinear operators q^'^q and their source, a back- 
ground U{1) baryon number gauge field [25-27]. We in- 
troduce a chemical potential through the U{1) baryon 
number gauge field which enters the DBI action in Ein- 
stein frame as 



ds^: 

where 



Kjr) 
E? 



K{r) 



_R2 

K{ry 



-dr^ 



i?2 



dxl + R^dnl , (6) Sd7 = -Tj I d^^e-^- det {e*/^P[GU + {2na')Fab) 



We allow a chemical potential through At{p) ^ 0. So 
the Action becomes 



rH = ttR^T. 



(7) 



S 



D7 



Witten identified this as the thermal description of the 
gauge theory in [4]. The parameter is of dimension 
one in the field theory and preserves the SO (6) symmetry 
so is identified as shown with temperature T. The black 



d X dp C 



-T-, 



d^x dp p{p)p\ 9t9l{l + L'^) 



9l 



A't^. 
(13) 
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where = (27rQ;')A(. In our convention of the metric 
this is 



(14) 



Now we can Legendre transform the action as we have a 
conserved quantity, the density, d ^fff-^ ■ 



where 



(15) 



We can redefine d = {27ra')T7d to give the simpler 
expression 



C = -Ty- 



f 4 4 
(W — W 



(16) 



By varying the Lagrangian with respect to A'^ , we get an 
expression for d{A[) which we can invert for an expression 
for ij(J) 



A' = d- 



(17) 



This can be used to find the chemical potential p ~ (2Tra') 



~ [W — Wjj I 



(l+L'2) 



\ /3(P) (»*+»!,) ^ P \ J 

(18) 

where p{p pn) — 0. 

The free energy can be found by integrating the Legen- 
dre transformed Lagrangian, the grand potential by inte- 
grating the original Lagrangian, where we replace A'^ (d) . 



5* 



157 



dp 



(u;4) 



^/3(p)Vl + i'2 



d2 



p^ {W^+Wjj) 



The grand potential is 



^ 7 JpH 



4 1 4 \ 2 



(H-i'2) 



where we need to note that F{p — >■ oo) = Q,{p — ?> oo) = 
p^, so we need to subtract ■|(Aj/y)^ from both integrals 
to renormalize them. 



III. ANALYSIS AND RESULTS 

The methodology to study the phase diagram of our 
model is straight-forward if laborious. We will work 
throughout in the massless quark limit. We can think 
of the scale A in the dilaton ansatz as our intrinsic scale 
of the theory and so we will leave that fixed. Then for 
each choice of parameters in the dilaton profile {A, T) we 
analyze the theory on a grid in T and p space. 

For each point on the T, p grid we seek three sorts of 
embedding. The flat embedding L = exists in all cases 
and describes the theory with jti = and c = 0. We use 
(18) to compute the d-p relation for these embeddings. 

We can also seek curved embeddings that miss the 
black hole. These solutions must have d = but are 
consistent for any value of p. Here we use the equation 
of motion for L from (16) and numerically shoot from 
an initial condition at p = with vanishing p derivative, 
L'(0) = 0. We seek solutions that approach L = at 
large p. These configurations have a non-zero conden- 
sate parameter c. 

Finally we can look for solutions that end on the black 
hole horizon. To find these we fix the density d and shoot 
out from all points along the horizon seeking a solution 
that approaches L = at large p. We then use (18) to 
compute p from the solution. In this way we can fill out 
the T,p grid. The condensate can again be extracted 
from the large p asymptotics of the embedding. 

After finding as many such solutions as exist at each 
point the easiest method to identify the transition points 
is to plot the density against p on fixed T lines. The tran- 
sitions and their order are then manifest. We display four 
sample plots in Fig 5 taken from scenarios below show- 
ing the four cases of the chiral transition and the meson 
melting transition being respectively first or second order 
in all combinations. 



A. Dependence on the change in coupling 

Let us first consider how the phase diagram depends 
on the height of the step in the gauge coupling function 
jS. We fix A (the intrinsic scale of the theory) and also 
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(a)r = 0.4, A = 30, r = 0.1, A = 1.715 
- shows a single first order transition. 




(b)T = 0.3, A = 3,r = 1, A = 
1.715, a = 3 - shows a first order 
transition followed by a second order 
transition as /i increases. 




(c)T = 1.55, A = 30, r = 1, A = 1.715 

- shows a second order transition 
followed by a first order transition as 
II increases. 




(d)T = 0.7, A = 3,r = 1, A = 
1.715,5 = 3 - shows two second order 
transitions. 

FIG. 5: Plots of density d versus chemical potential /i. The 
top lighter line (green) in each case corresponds to the flat em- 
bedding; the horizontal line (black) along the axis is a chiral 
symmetry breaking (Minkowski) embedding; the near vertical 
dark line (blue) is a black hole embedding. Transition points 
are shown by the dotted vertical lines. 



r = 1 and explore the phase structure as a function of 
A. We display the results for three choices of A in Fig 6. 

In these and all our future phase diagrams the regions 
shown are similar to those in Fig 2(a) we will simply 




0.5 1.0 1.5 2.0 2.5 3.0 3.5 



(ii)A = 3 




ic)A = 30 

FIG. 6: Plots for three possible phase diagrams for the 
choices ^ = 3, 15, 30. Large (small) A gives second (first) 
order transition at low T. F = 1, A = 1.715. 



display the phase boundaries and their order henceforth. 

As mentioned in footnote 2 above there is a critical 
value of A for chiral symmetry breaking to occur. A 
conformal theory can not break a symmetry since it of- 
fers no scale for that symmetry breaking to occur at. In 
fact some finite departure from conformality is needed 
to break the chiral symmetry. For these choices of A, F 
the critical A is Ac = 2.1. We work above this value 
throughout. 

At low A there is a single transition for chiral symmetry 
restoration and meson melting which is first order for all 
T and /i. On the gravity side this is a transition between 
the curved embedding that misses the black hole and the 
flat embedding. In this case an embedding ending on the 
black hole never plays a role. 

For larger A, a new phase with chiral symmetry break- 
ing but melted mesons develops. There is a regime now 
in which the curved embedding ending on the black hole 
is energetically favoured. The transition from the chiral 
symmetry breaking phase to this new phase is second or- 
der. The chiral symmetry restoration phase remains first 
order. 

At very large A the chiral restoration transition be- 
comes second order at high density. This latter phase 
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{a)T = 0.1 



2 3 4 5 

(b)T= 1.0 



action cost of entering the region in which the dilaton 
is large. As temperature is introduced through a smah 
horizon the interior of the space is "eaten" but the D7 
embedding remains obhvious to this change since it never 
reaches down to small r. As temperature rises the point 
of transition is when the horizon moves through the scale 
A where the dilaton step is. Once the region with a large 
dilaton is eaten by the black hole the preferred D7 em- 
bedding is the flat one. 

In Fig 7 we show an extreme case of this behaviour 
explicitly. Here we have taken T very large so that the 
transition in the dilaton between the low and high value 
is very sharp. We plot the f3 profile against our radial 
parameter w and mark in red the position of the black 
hole horizon. Note that in the w coordinates the region 
where /3 is large depends on the temperature (it doesn't 
in the orignal r coordinate) . The dilaton effective radius 
A, is 



A* — 



'A2 + VA4 - 4r4 



(20) 



(c)r = 1.20208 = Dilaton is 
almost screened. The residual effect is 
due to finite F effect. 



FIG. 7: Plots for parameter choices A = 5,r = 100, A = 1.7. 
The blue lines show the value of the coupling /3. The red line 
shows the position of the horizon. The final plot corresponds 
to the point of the first order transition. 



resembles that of the theory with chiral symmetry break- 
ing induced by a magnetic field [15]. In fact the B field 
case can be thought of as our case but with a choice of /3 
given by 




(19) 



It is the black dotted curve {wh = 0) in Fig 4 - it is not 
surprising therefore that we see similar phase structure 
here (and indeed that we do provides strength to our 
analysis which is capturing the behaviour of top down 
models). 

For very large A the step becomes very sharp and there 
is little change relative to our phase diagram in Fig 6. 
In particular the thermal transition always remains first 
order. 

The behaviour we are seeing here can be readily ex- 
plained from the D7 perspective. First of all the zero 
density transition with temperature is first order for a 
simple reason. The D7 embedding breaks chiral symme- 
try at zero temperature because it prefers to avoid the 



where the argument of tanh in (12) vanishes. So, as T in- 
creases A* decreases. When T becomes Tc = A, = Tc 
the dilaton is perfectly screened by the black hole hori- 
zon, (i.e. If T = 0, A* = A. If T = ^, A* = ^). The 
point of the first order transition is where the horizon 
screens the dilaton. 

When density is introduced the story can become more 
complex. The action is (16) where it can be seen from the 
first of the two terms in the square root that including 
d increases the action. This increase can be beneficial 
though if the second term with /3 can be reduced. It is 
possible to reduce the /3 term if the D7 enters the region 
where /3 is large at small p. This means that the situation 
can arise where curving off the axis and then spiking on to 
the axis can be the lowest action state. This is typically 
more likely where 13 is largest in the interior space and 
the most savings can be made entering that region at low 
p. As we have seen at large values of A cmbcddings that 
spike onto the horizon do play a role introducing an extra 
phase. 

It is only possible to have second order transitions if 
all three phases we have described are present. In the D7 
description the D7 must move from a curved embedding 
that avoids the black hole to a configuration that spikes 
onto the black hole to a flat embedding smoothly. 



B. Dependence on the speed of running 

The parameter F controls the period in p or RG scale 
over which the change in the coupling A occurs. It allows 
us to naively go from a precociously running theory to a 
walking theory (although the change in the parameter A 
over that period may enter into what is meant by walking 
versus running too). 
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FIG. 8: Example plots of three possible phase structures for 
A = 30, A = 1.715 and varying F. Large (small) F gives a 
second (first) order transition at low T. 



In Fig 8 we show the phase diagram as a function of 
r at fixed A and A. We start at F = 1 with a configura- 
tion already discussed that has all three phases present 
and second order transitions at high density. As F is re- 
duced so that the step function in the dilaton becomes 
broader the first order nature of the transitions reasserts 
itself. By F = 0.1 the mixed phase with chiral symmetry 
breaking but melted mesons is no longer preferred at any 
temperature or chemical potential value - there is a single 
first order transition. 

In conclusion then moving towards a walking theory by 
either increasing the width of the running or decreasing 
the magnitude of the increase in the coupling both move 
us towards a first order chiral transition. Stronger or 
quicker running favours a second order transition at low 
temperature, high density. 



IV. BREAKING THE p-L SYMMETRY 

Our goal is to attempt to reproduce a phase diagram 
comparable to that of QCD in our holographic model. So 
far we have failed to generate a second order transition 



FIG. 9: Sample phase diagrams for theories with none zero 
a. 



with temperature at zero density which is a key part of 
the QCD picture. 

We have a further natural freedom within our holo- 
graphic model to exploit though. Our running dilaton is 
in someway supposed to represent the backreaction of the 
quark fields on the strongly interacting gauge dynamics 
to allow us to model theories with more interesting dy- 
namics than the conformal A/" = 4 gauge fields. We in- 
troduce quarks through D7 branes that break the 5*0(6) 
symmetry of the five sphere of the original AdS/CFT 
Correspondence down to 50(4) x 50(2). Our metric 
and ansatz for the running coupling (2) though respected 
the full 5*0(6) symmetry. It seems reasonable to make 
use of the broken symmetry to introduce a further free 
parameter into our model. 

The most succesful scenario we have found is to intro- 
duce our explicit L — p symmetry breaking parameter, a 
through the blackening factors of the metric 



gt 



(«;4 



w 



4 \2 

h) 



W 



H 



with 



— ?► p'^ H , a > 1. 

a 



(21) 



(22) 
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We show the a dependence of our model in Fig 9. We 
start from a model with a first order transition through- 
out the phase plane. As we increase a a region with 
melted mesons but chiral symmetry breaking develops, 
associated with second order transitions. It then spreads 
to the zero chemical potential axis 

The reason for the onset of second order transitions 
with just temperature is simply understood. We have 
deformed the black hole horizon into an ellipse whose 
major axis is along the L axis. The area of enlarged 
dilaton remains circular in the p-L plane. Thus there are 
temperature periods in which the area of the p-L plane 
with a large dilaton is covered except for a small piece 
that emerges from the horizon near the p axis. If the 
value of the dilaton is sufficiently large in that uncovered 
area to encourage the D7 to avoid it, but the horizon on 
the L axis has met the zero temperature D7 embedding, 
then a second order transition to a black hole embedding 
is likely. Since in the absence of the rise in the dilaton the 
flat embedding would now be preferred the D7 settles on 
the horizon so it just misses the raised dilaton area. As 
the black hole grows further the embedding is likely to 
track down onto the axis smoothly as the raised dilaton 
area is finally eaten. This intuition is indeed matched by 
the solutions as shown in Fig 9. 

The bottom phase diagram in Fig 9 achieves our goal 
of reproducing a chiral transition that is second order 
with temperature but first order with density. 

V. SUMMARY 

In this paper we have converted the D3 probe-D7 sys- 
tem, that holographically describes A/" = 4 super Yang- 
Mills theory with quenched N = 2 quark multiplets, to a 
phenomenological description of strongly coupled quark 
matter. Wc introduced a simple unback-reacted profile 
for the dilaton that describes a step of variable height 
and width in the running coupling of the gauge theory - 
(2). This breaks the conformal symmetry of the model 
and introduces chiral symmetry breaking. We have then 
studied the temperature and chemical potential phase 
structure of the model. 

The phase diagrams consist of three phases - a chi- 
rally symmetric phase at large temperature and density; 
a chirally broken phase with non-zero quark density at 
intermediate values of T and p; and a chiral symmetry 



broken phase with zero quark density at low T and p. Fig 
2(a) shows these phase and their holographic analogue in 
a previously studied case where chiral symmetry was bro- 
ken by an applied magnetic field. Here we showed that a 
small wide step in the gauge coupling's running gives rise 
to a single first order transition between the chiral sym- 
metric and the broken phase (see Fig 6). If the step is 
made larger in height or thinner then the chirally broken 
phase with non-zero density also plays a role. Here the 
transitions at low temperature with chemical potential 
can be second order. These results match known results 
in top down models in the presence of magnetic fields to 
induce the symmetry breaking. 

We were interested in reproducing phase diagrams with 
the structure believed to exist in QCD. To do this we 
made use of the broken 5*0(6) symmetry of the gravity 
dual in the presence of D7 branes. Were the branes back- 
reacted the dilaton and geometry would reflect this sym- 
metry breaking. We introduced a further phenomenologi- 
cal parameters a in the black hole blackening factor. This 
models allowed us to control which volumes of the holo- 
graphic space have a large dilaton value within, which the 
D7 branes prefer to avoid. Using this one extra parame- 
ter we were able to generate phase diagrams like those in 
QCD with a chiral restoration transition that was second 
order with temperature but first order with density (see 
Fig 9). 

The ease with which such a variety of phase struc- 
tures could be obtained is very encouraging for the idea 
of phenomenologically modeling the QCD phase diagram 
holographically. Further the phenomenological parame- 
ters we introduced are very natural in this context and 
it seems likely that top down models with such phase 
structures should be possible as a result. 
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